The Radius of the Proton: Size Does Matter 



Jonathan D. CarrolP, Anthony W. Thomas^, Johann Rafelski'' and 

Gerald A. Miller^ 

'^Centre for the Subatomic Structure of Matter ( CSSM), 
Department of Physics, University of Adelaide, SA 5005, Australia, 
http://www.physics.adelaide.edu.au/cssm 
^Departments of Physics and Mathematics, University of Arizona, Tucson, Arizona, 85721 USA 
"University of Washington, Seattle, WA 98195-1560 USA 

Abstract. The measurement by Pohl et al. [1] of the 2S^^' to 2P'^j^ transition in muonic hydrogen 
and the subsequent analysis has led to a conclusion that the rms charge radius of the proton differs 
from the accepted (CODATA [2]) value by approximately 4%, leading to a 4.9 a discrepancy. We 
investigate the muonic hydrogen spectrum relevant to this transition using bound-state QED with 
Dirac wave-functions and comment on the extent to which the perturbation-theory analysis which 
leads to the above conclusion can be confirmed. 

Keywords: proton radius, QED, lamb shift 



INTRODUCTION 

In this work we calculate the transition energy relevant to the aforementioned experiment 
of Pohl etal. [1] (as depicted in Fig. 1) using the Dirac equation in an attempt to quantify 
the errors associated with the perturbative approach. In the sections following, we 
discuss the nature of the transition and its components; the method by which we calculate 
the energies corresponding to the various eigenstates; and the predicted energies of the 
component shifts as a brief account of a longer upcoming publication [3] in which we 
shall detail the components in full with comparisons to previous work [4, 5, 6]. We note 
that since this talk was presented, we have investigated an additional term that to our 
knowledge does not already appear in the analysis of Pohl et al. and which may account 
for all or part of the discrepancy [7]. 

NUMERICAL METHOD 

To calculate the theoretical shift corresponding to the measured transition, previous 
authors have primarily used perturbation theory with non-relativistic wave-functions 
to predict the size of the contributing effects, including relativistic effects. To better 
approximate the exact energies, we can use the Dirac equation for the muon with the 
appropriate potential as an effective approximation to the two-particle Bethe-Saltpeter 
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FIGURE 1. (Color online) Muonic hydrogen spectrum, including finite-size correction. Lamb shift, fine 
structure, and hyperfine structure. Also shown is the measured 28^^"^' to ^V^J^ transition (green, dotted, 
^) as per Ref. [1]. 



equation [8] to calculate the perturbed wave-functions, expressed here as a spinor 
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normalised to unity, such that the probability is 

r^[ga{rf + fa{rf] dr 



noting that 
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Since this is a relativistic system, we use the reduced mass ji in place of the muon mass 
in the Dirac equation 
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which along with the addition of recoil corrections provides a good approximation to 
the Bethe-Saltpeter equation. Since the binding of the muon in this system is extremely 
weak, the eigenvalue Ca for each state calculated using the Dirac equation is approxi- 
mately equal to the reduced mass ji. In order to precisely calculate the variance from 
this value, we shift our eigenvalue down by the reduced mass, such that the eigenvalue 



we are now solving for is = — /i, thus the effective Dirac equation is 




where the value of Ka is specific to each eigenstate, namely 

1Si/2:Jc= -1, 2Si/2:Jc= -1, 
2Pi/2:Jc= +1, 2P3/2:Jc= -2. 

The (shifted) eigenvalues can be reliably reproduced by using the point-Coulomb 
potential 

V{r) = - — , (6) 
r 

in Eq. (5). In order to integrate Eq. (5), we supply an initial guess for the eigenvalue 
Xa, and appropriate boundary behaviour of upper and lower components of the wave- 
function at small and large radii, then integrate from each limit towards a central match- 
point. The discontinuity in the wave-function integrated from each limit is used as a 
measure of the inaccuracy of the eigenvalue and a refined estimate is calculated. This 
process is iterated until the change in Xa is less than the required tolerance, at which 
point we regard the wave-function to be converged. 

To convince ourselves that our method is self-consistently accurate, we check the 
accuracy of our procedure using several methods. The unperturbed Dirac eigenvalues 
are known analytically [9] to be 
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where is the principle quantum number for the state a. We first ensure that we are able 
to reproduce these values. For the 2S j /2 wave-function, we reproduce this value to within 
0.01 jUeV using quad-precision Fortran, a sufficiently large grid size, and sufficiently 
small grid spacing, within reasonable compute-time. We also check the validity of the 
virial theorem for our solutions (refer to Ref. [9] for further details) by calculating the 
reduced eigenvalue as 

X = (2Si/2 I juj8 +V{r)+r- W (r) 1 2S1/2) - jU, (8) 

— * 

which tests the accuracy of the wave-function at the origin where |Vy| is greatest. We 
calculate that the values obtained using Eqs. (7) and (8) differ by 0.18 jieV for a point- 
Coulomb potential, and 0.45 jUeV a finite-Coulomb potential (to be discussed later). We 
therefore conservatively take our errors to be of the order of ~ ±0.5 jUeV. 



NUMERICAL CALCULATIONS 



• 2Si/2~2Pi/2 Lamb Shift: The Lamb shift is the splitting of the otherwise degener- 
ate 2Si/2 and 2Pi/2 eigenstates attributed to the vacuum polarization potential Vyp 
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We can calculate the shift in eigenvalues using converged Dirac wave-functions in 
response to the Coulomb and vacuum polarization potentials, and in this case we simply 
take the difference between the converged eigenvalues for the 2Si/2 and 2Pi/2 states 

^S'' = ^2Pi/2-^2Si/2 = 205.1822meV. (10) 

Care must be taken when comparing this calculation to that of perturbative results since 
our calculation includes relativistic corrections, which are included later as corrections 
in perturbative calculations, e.g. Ref. [10]. 

• Proton Finite-Size Corrections: To calculate this effect in our fully relativistic 
calculation, we consider the replacement of the point-Coulomb potential with the finite- 
size Coulomb potential in Eq. (5) 

Vc{r) = -—^-Zajp^d\ (11) 
r J \r — r \ 

where p(r) is the proton charge-distribution (or more accurately, the slope of the electric 
form-factor). We have studied the dependence of the finite-size correction on the form 
of this term (always normalised to unity) and this will be summarized in an upcoming 
publication (Ref. [11]), though the dependence on the choice of charge-distribution — 
whether it be exponential, Yukawa, or Gaussian in form — appears to be small. Similarly 
the finite vacuum polarization potential is given by a convolution of Eq. (9) with the 
charge-distribution. 

The exponential form for the charge-distribution, normalised to unity such that 
/ p (r) d^r = 1 is given by 

^^"^^1^"^""' ^ = ^^V{rl). (12) 

We calculate the Lamb shift by taking the difference between the appropriate eigen- 
values calculated using the Dirac equation with the potential given by Eq. (11) with the 
charge-distribution given by Eq. (12) for various values of (r^). We then interpolate the 
energy shifts and fit the data to a cubic of the form 

f{x)^A + B{rl)+C{?/l\ (13) 

which provides the relevant parameterization. The dependence of the Lamb shift on the 
rms charge radius in the presence of an exponential finite-sized Coulomb potential and 



finite vacuum polarization potential is given by 



AEflnite = 205. 1822 -5.2519(?^) + 0.0546 (rj)^/^ meV. (14) 

• 2P Fine Structure: Subtracting the converged eigenvalues of the and 2P3/2 
eigenstates gives the fine structure splitting 

which we can also calculate in the presence of the various potentials. For the case of an 
exponential finite-Coulomb potential with finite vacuum polarization potential, the 2P 
fine structure shift is 

b£f^ = 8.4206(5) meV. (16) 

The shifts due to finite-size effects (as compared to the point cases) are below the 
level of errors for our calculation. The point vacuum polarization itself increases the fine 
structure shift by 5 jUeV. 

• 2S1/2 Hyperfine Structure: The splitting between the 28 F = and F = I 

hyperfine eigenvalues is given [12] by 

A£Sr'-"^"'=^/JrlV(0)p. (17) 

The value of the 2S hyperfine shift, as calculated using Eq. (17) with the wave- 
function calculated with the Dirac equation in the presence of the combined point- 
Coulomb and point vacuum polarization potentials is 

AEjj^P^ = 22.8967(5) meV. (18) 

The finite-size effects will be investigated in an upcoming publication [3]. 

• 2Pi/2 Hyperfine Structure: The 2Pi/2 Hyperfine structure is of no consequence to 
the measured transition we are investigating. Nonetheless, we calculate the energy of the 
2P^^2° and '2P\j^ states as a confirmation of our method, and to compare to perturbative 
results. The 2P hyperfine structure is given [12] by 

4fs = 2i8 (^^) {Fmp \ I • J | Fm^) , (19) 

where the non-zero terms in the dot-product are given by 

{Fmp 1 1 • J I Fmp) = \ [F{F + 1) - (20) 
which, for Schrodinger wave-functions gives 

^Zf = l^r/al , (21) 



to which anomalous magnetic moments provide further corrections. Using the converged 
Dirac wave-functions with exponential finite-Coulomb and finite vacuum polarization 
potentials (rather than Schrodinger wave-functions) we calculate the expectation value 
of and find 

AEjjpg/' = 7.6204(5) meV. (22) 

The addition of the (point) vacuum polarization potential to the point-Coulomb potential 
increases the splitting by 0.0017(5) meV, and the introduction of the finite-Coulomb po- 
tential increases this further by 0.0045(5) meV. The finite vacuum polarization potential 
does not alter the result from the point case here. 

• 2P3/2 Hyperfine Structure: Following the same method as the 2Pi/2 calculation, 
we can calculate the energy levels for the ^"^^ -^^s/i^ eigenstates. Using the 

converged Dirac wave-functions we find 

AEjJg/' = 3.0415(5) meV (23) 

when the potential consists of the exponential finite-Coulomb and finite vacuum polar- 
ization potentials. For this state, the addition of the (point) vacuum polarization potential 
to the point-Coulomb potential increases the splitting by 0.0007(5) meV, and the intro- 
duction of the finite-size effects was found to make no change within the limits of our 
calculation. 

CONCLUSIONS 

We find that the Dirac calculations performed here agree well with perturbative results 
once appropriate corrections are made (taking care regarding double-counting of ef- 
fects). The calculations presented here and discussions of the comparison to perturbative 
calculations will be fully detailed in several upcoming publications [3, 11]. 
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